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1.  INTRODUCTION 


Conventional  classical  supercomputers  will  likely  reach  their  computational  and/or  feasible 
commercial  cost  limit  in  20-30  years.  It  is  extremely  unlikely  that  many  of  the  most  difficult 
and  important  non-linear  field  theories  like  Navier-Stokes,  QCD,  General  Relativity,  and 
quantum  gravity  will  be  solved  before  these  limits  are  realized.  By  far  the  most  promising 
way  forward  appears  to  be  through  quantum  computer  development.  It  is  also  generally  ac¬ 
knowledged  that  quantum  computers  will  provide  the  most  powerful  system  for  simulations 
of  quantum  mechanics. 

This  paper  develops  primarily  type-II  quantum  computer  lattice-gas  algorithms  because  of 
the  apparently  formidable  unsolved  technological  obstacle  of  establishing  completely  globally 
phase-coherent  type-I  quantum  computation.  In  particular  the  Weyl  ID,  2D,  and  3D  and 
Dirac  ID  lattice-gas  equations  for  relativistic  quantum  point  particles  are  given. 

Section  two  contains  a  derivation  of  the  orthonormal  properties  of  the  W  matrices.  Section 
three  discusses  the  lattice-gas  algorithms  of  the  ID  Weyl  and  Dirac  equations  in  the  context 
of  the  papers  of  [Bialynicki-Birula,  I.,  1994]  and  [Meyer,  D.  A.,  1996].  Explicit  lattice-gas 
algorithms,  first  order  solutions,  and  implementations  are  presented  for  a  specific  represen¬ 
tation.  In  section  four  an  explicit  first  order  lattice-gas  solution  of  the  2D  Weyl  equation  in 
a  specific  representation  is  presented.  The  exact  two  component  spinor  solution  to  the  2D 
Dirac  equation  is  given,  A  numerically  untestedresult  and  possible  methods  for  obtaining 
the  2D  Dirac  lattice-gas  solution  is  included.  Insection  five  an  explicit  first  order  lattice-gas 
solution  to  the  3D  Weyl  equation  is  presented  in  a  specific  representation. 
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2.  ORTHOGONALITY  AND  UNITARITY  CONDITIONS  OF  THE  DISCRETE 
WEYL  AND  DIRAC  PROPAGATOR 


The  Weyl  equation  describes  the  spinor  of  a  rest-massless  relativistic  quantum  point  object 
while  the  Dirac  equation  describes  the  spinor  of  a  massive  relativistic  quantum  point  object. 
For  the  Weyl  solution,  [Bialynicki-Birula,  1994]  postulates  a  two  component  spinor  <f>  on 
a  manifold  with  a  body-centered  3D  discrete  Euclidean  space  lattice  and  a  Newtonian  time 
with  the  following  Huygens-like  propagation  for  each  discrete  update  time  r.  W(h)  is  a  2x2 
matrix  operating  at  space  point  r  +  h, 

<f>(r,t  +  r)  =  W(h)<f>(r  +  K,t ),  (1) 

it 

but  for  the  Dirac  solution  W  becomes  a  4x4  matrix  D  and  (f>  becomes  a  four  component 
spinor  <fi  1  .  Each  of  the  members  of  the  set  of  eight  primitive  lattice  vectors  h  start  at  f 
and  end  at  one  of  the  eight  corners  of  the  cube  defined  by  the  first  octant  of  the  coordinate 
system  originating  at  r. 

To  determine  the  algebra  of  the  2x2  W  matrices  first  consider  two  equal-time  independent 
spinors: 

<Kfut  +  r)  =  J2W(hiWl  +  h,t)  (2) 

hi 

^(ra.*+V)  =  Y.W\h2)4>\f2  +  h2,t). 
h2 

Now  consider  the  overlap  sum: 

£  <t>](f2,  t  +  r)4>{f\,t  +  r)£r- +/lVi+£  (3) 

ri  ,r1 

Thus  for  set  elements  satisfying  h  ^  h » the  above  Kronecker  delta  function  allows  the  prod¬ 
uct  of  partially  overlapping  (  r\  ^  r2  )  spinors  to  be  non-trivial.  But  we  will  demand,  as 

XA11  tilded  operators  following  are  4x4  matrices  and  non-tilded  matrices  are  2x2. 
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Bialynicki-Birula  did,  that  the  sum  of  these  products  is  trivial.  The  Kronecker  delta  also 
defines  a  maximum  separation  for  non- trivial  overlap  that  occurs  when  h  =  -ii\  e.g.  at 
opposite  comers  of  the  first  octant.  Equivalently,  and  using  Eq.  (2): 


Hence; 


?2,  i  +  r)<p(fi,  t  +  h'  = 

r'lj  T2 


E  Wt{h2)W(h1)(t>Hf2  +  h2.,t)<f)(f1  4-  hi,t)6r-.2rl+K^ 


E^^i  +  h-  h',t  +  T)(f)(fi,t  +  r)  = 

Y  ^\h2)W{hi)(t>\f{  +  h-ti  +  h2,t)4>(fi  +  hut). 

The  normalization  condition, 


Y  4>Hr,t)(f>(f,t)  =  1 

T 

can  be  re-written  for  the  new  time  t  4-  r  as 

Y  +h-ti,t  +  r)^(rj,t  +  r)  = 

A 

and  hence: 


(4) 

(5) 

(6) 

(7) 


Y  ^f(rl  +  h-h'  +  h2, t)<f>(fi  +  hut)  =  6n_ 


h'+h2,hi. 


(8) 


Thus  ,  we  have  equivalently  from  Eqs.  (5),  (7),  and  (8)  the  orthonormaJ  conditions  for  the 
W  matrices:  « 5_f,  0  =  Z  (9) 

huh? 


And  equivalently  the  orthonormal  conditions  are  more  simply: 

k-h>,o  =  E  W'(hi  +  h'~  h)W(hi). 

hi 


(10) 
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3.  THE  ONE-DIMENSIONAL  WEYL  AND  DIRAC  EQUATIONS 


Bialynicki-Birula’s  ID  Weyl  algorithm  and  Meyer’s  ID  Dirac  algorithm  [Meyer,  D.  A.,  1996], 
e.g.  Eq.  (1),  both  have  the  same  implicit  form, 

<t>(x,t  +  T )  -  W(-h)^)(x  -  h,t)  +  W{h)<j)(x  +  h,t)  (11) 

but  the  explicit  matrix  elements  are  different.  Though  it  is  not  by  design,  Meyer’s  matri¬ 
ces  (equation  (17)  of  [Meyer,  D.  A.,  1996])  substituted  for  Bialynicki-Birula’s  ID  matrices 
reproduces  the  required  orthonormality  and  Bialynicki-Birula’s  helicity  relation, 


W  (~ll)  —  (Ty  W*(h)  Oy. 


(12) 


The  Bialynicki-Birula  ID  algorithm  (see  Eqs.  (39)  and  (40)  of  [Bialynicki-Birula,  I.,  1994]), 
for  the  Dirac  equation  yields  the  same  correct  result  as  Meyer’s  because  the  four-spinor 
difference  update  equations  decouple  into  two  independent  particle  and  anti-particle  two- 


component  spinor  update  equations,  each  equivalent  to  Meyer’s  result.  Writing  6  =  {  f  L 

\<PR 

expanding  and  collecting  components  the  difference  update  equations  are: 


/  <£x,(x,t  +  r)\  /  Cos(0)  iSin(6)  \  (  <f>  L(x  +  h,t)\ 

\  0  t  +  T)j  \  iSin(9)  Cos{6)  )  \  <t>  r(x  -  h,t) ) 

Now  by  using  a  streaming  operator  Sx, 


a  (<f>L(x,t)\  =  (<t>L(x  +  h,t)\ 

x\4>R(x,t)J  \  <t>  R(x  -  h,t)  ) 

one  can  re-write  the  lattice-gas  equation  (13)  as: 


(14) 


/  (, pL(x,t  +  T ) 

\  <t>  R(x,t  +  T ) 


with  C  given  by  the  collision  operator  shown  in(13)  above.  UsingMathematica]  Wolfram, S., 
2001]  this  lattice-gas  algorithm  can  be  implemented  using  symbolic  terms  which  can  then 
be  expanded  in  terms  of  the  lattice  spacing.  The  following  symbolic  psuedo-code  difference 
scheme  has  lattice  period  r  =  2  and  the  amplitudes  <f>  are  initialized  at  time  t=0: 


if  mod[n,2]  =0; 

<j>L[l,n]  =  Cos(0)(f>L[l,n  —  1]  +  iSin(0)(j>R[l,n  —  1]; 
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0 i?[^>  w]  iSin[6^cpi,[l ,  n  1]  +  Cos(6^cpji[l,  n  —  l]j 

if  mod[n,2]  =  1; 

<t>L[l,n]  =  4>l[1  +  1  ,n  —  1]; 

<j>n[l,n}  =  <pR[l-l,n-l); 
do  l,l,numberofgridpoints 
<Pl[1,  0]  -phiL[l]- 
0fl[^,O]  =  phiR[l}; 
end  do 

Once  the  above  has  been  executed,  the  Mathematica  [Wolfram,  S.,  2001]  command 

Table[{(t>L[l,  n],<f>R[l,  n]},  {t,  0, 2}],  (16) 

or  some  equivalent  form  of  Do  loop  then  produces  the  values  of  <p  for  all  the  integer  t  values. 
These  terms  then  must  be  expanded  to  first  order  in  the  lattice  spacing  using  fourier  series 
(see  equation  (24)  for  an  example). 

For  6  =  0  Bialynieki-Birula’s  Weyl  equation  algorithm  is  recovered  and  to  first  order  in 
the  lattice  spacing  <5x  (e.g.  use  approximations  like  <fL[l  - 1, 0]  =  <f>L[l,  0]  4-  8x<j>{l’0)[l,  0]),  this 
reproduces  the  ID  massless  Weyl  equation.  Bialynicki-Birula’s  ID  Weyl  matrices  can  be 
explicitly  found  from  linear  combinations  of  his  P  matrices:  W(h)  =  (Pi  -f-  P3)/2  and  W(-h) 
=  {Pi  -  -P4)/2.  For  the  case  6  =  me  where  e  is  an  infinitesimal  scalar  the  ID  massive  Dirac 
equation  results  to  first  order  in  both  e  and  the  lattice  spacing.  This  is,  of  course,  what  Feyn¬ 
man  proved  [ Feynm.an ,  R.  P.  and  Hibbs,  A.  R.,  1965].  And,  as  required,  the  Weyl  equation 
results  for  m=0  and  this  is  mathematically  equivalent  to  the  case  6  =  0.  Although  there  ap¬ 
peared  no  recognizable  symbolic  result  for  a  finite  and  non-trivial  6  Meyer  gives  some  results 
[Meyer,  D.  A.,  1996]  connecting  non-trivial  finite  6  values  with  finite  velocity  propagation. 


The  Weyl  equation  and  the  Dirac  equation  can  thus  be  modeled  with  the  same  algorithm 
but  the  Weyl  model  requires  all  finite  matrix  parameter  values  and  the  Dirac  model  requires 
some  infinitesimal  matrix  parameter  values.  Bialynicki-Birula’s  stated  algorithms  for  the 
ID  Weyl  and  Dirac  equations  are  correct.  Since  it  is  not  exactly  clear  how  to  interpret  the 
Dirac  algorithm  using  non-infinitesimal  6  values  one  must  choose  some  e  magnitude  that 
is  sufficiently  small  for  the  precision  required  to  reproduce  the  behavior  of  a  Dirac  object. 
Bialynicki-Birula’s  2D  and  3D  algorithms  behave  similarly. 

One  can  derive  the  ID  and  2D  Dirac  differential  equations  and  the  ID,  2D,  and  3D  Weyl 
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differential  equations  by  considering  Dirac’s  equation 

he 

ihd  =  — —ot  •  Vip-hmc2/3ip  (17) 

i 

where  a  and  (3  are  each  replaced  by  one  of  the  2x2  Pauli  matrices  o’  and  m=0  for  the  Weyl 
equation.  Thus  there  are  six  possible  representations,  or  equations,  for  the  ID  Dirac  equation 
and  three  possible  representations  for  the  ID  Weyl  equation.  The  'square'  of  each  equation 
yields  a  Klein-Gordon  equation.  Only  for  the  massive  case  in  three  space  dimensions  are  4x4 
matrices  required  for  equation  (17).  [Jacobson,  T.  and  Shulman,  L.  S.,  1984]  take  a  =  crz 
and  P  =  —ox  which  gives 

he 

ihd  tip  i= — r -dxipi — mc2ip2 

(18) 

Tic 

ihdtip 2  =  — r-  dtip 2 — mc2ipi 
% 


These  equations  are  solved  to  first  order  in  both  the  lattice  spacing  and  in  e  by  the  lattice-gas 
algorithm  (15)  and 

c=  (  cos(me)  isin(me)  \ 

y  isin(me)  cos(me)  J  ^  ' 
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4.  THE  TWO-DIMENSIONAL  WEYL  AND  DIRAC  EQUATIONS 

As  with  Bialynieki-Birula  we  choose  =  1/2  Pi  so  the  2D  W  matrices  are 


W,  = 


0) 

0  r 


W. 


0 

0 


w2  = 


w. 


(20) 


0  -1 
0  1 


The  2D  Weyl  differential  equation  is  obtained  from  Dirac’s  equation,  (17),  with  m=0: 

%Q 

ihdt$=  (a  xdxfa+  aydyip).  (21) 

There  are  six  possible  representations,  or  equations,  depending  on  which  components  of 
are  substituted  for  which  of  oix  and,  cty .  One  possible  choice  olx  —  <r2,  ocy  =  gives  the 
explicit  equation 


dti> i  =  c{~dxtp i  -  dyip2 ) 


dti) 2  =  c{+dxi)2  -  dyip i). 


(22) 


The  following  unitary  lattice-gas  algorithm  reproduces  this  representation  of  the  2D  Weyl 
equation  to  first  order  in  lattice  spacing: 


/  <t>L(x,y,t  +  r ) 
\4>n{x,y,t  +  T ) 


(0>  Sy  Ca 


q  \( (pL(x,y,t) 
**)[<{>  R(x,y,t) 


(23) 


The  explicit  result  after  one  period  r  is  given,  using  Mathematical  Wolfram,  S.,  2001],  by 


i(cpL[l,  m,  r)  -  4>L[l,  m,  0])/e  =  [l,  m,  0]  +  ^,0,0)  [l,  m,  0])/e,  (24) 

where  Sx  and  Sy  act  as  the  streaming  operator  S  from  before  but  now  act  only  on  the 
indicated  variable.  The  C  matrices  were  determined  algebraically  by  expanding  out  the 
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difference  equation  (21)  and  equating  the  corresponding  matrix  elements  of  relation  (1)  sub¬ 
stituted  with  the  above  2D  W  matrices  in  (20): 


=  'i1)“da=7i(-i  0- 

Then  the  helicity  relation  (12)  yields 


(25) 


W  {h) 

W  (-h) 

Wi 

w4 

w2 

w3 

w3 

w2 

w4 

Wi 

Substitution  of  these  2  x  2  W  matrices  into  equations  (39)  and  (40)  of  [Bialynicki-Bir-ula,  I. 
1994]  gives  a  four  component  difference  equation  for  the  2D  Dirac  equation. 

The  2D  Dirac  equation  is  obtained  from  Eq.  (17): 


he 

ihdt ij)  =  —(axdxip  +  otydyi/j)  +  /3mc2ip. 
% 


(26) 


There  are  3!  possible  representations,  or  equations,  depending  on  which  components  of  a*  are 
substituted  for  which  of  as,  ay  and  (3.  One  possible  choice  ax  —  crx,  av  —  oy,  and  /3  =  a z 
gives  the  explicit  equation 


Tic 

ihdtipi  =  —(dxtj>2  -  idyipz)  -f  mcVi 
% 

Tic 

ihdti>2  =  —(dxipi  +  idyl/)!)  -  mc2V;2- 
% 


(27) 


The  formal  solution  to  Eq.  (26)  is  obtained  by  determination  of  S  for 

%b'(x  ,y,t')  =  Si/>(x,y, t), 


(28) 
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with  the  result: 


5  = 


E  +  me 2 
2  me 


/  1  _£=£ . 

I  £4-mc2 

1 

£+mc2  1 


(29) 


Bialynieki-Birula’s  4x4  D  matrices  for 


<P(^  t  +  r)  =  23  D(h)4>(r  +  h,  t) 

h 


are  explicitly  represented  by: 


Cos  (me) 

0 

0 

—iSin(me)  \ 

(  0 

Cos(me) 

iSin(me ) 

o\ 

1 

Cos(me) 

0 

0 

iSin(me) 

£  1 

0 

Cos(me) 

—iSin(me) 

0 

2 

iSin(me) 

0 

0 

—Cos(me) 

’V2~  2 

0 

iSin(me) 

Cos(me) 

0 

iSin(me) 

0 

0 

Cos  (me)  J 

u 

iSin(me) 

—Cos(me) 

0/ 

Cos(me) 

0 

0 

iSin(me )  \ 

0 

—Cos(me) 

iSin(me) 

0\ 

1 

—Cos(me) 

0 

0 

iSin(me ) 

Dd-i 

0 

Cos(me) 

iSin(me) 

0 

2 

iSin(me) 

0 

0 

Cos  (me) 

’  Ua  ~  2 

0 

—iSin(me) 

Cos(me) 

0 

\-iSin(me) 

0 

0 

Cos(me)  / 

VO 

iSin(me) 

Cos(me) 

0/ 

(31) 


D, 


Since  Bialynicki-Birula’s  difference  equation  for  the  2D  Dirac  equation  uses  four-component 
spinors  but  the  solution  of  the  2D  Dirac  equation  mathematically  only  necessitates  two- 
component  spinors,  a  lattice-gas  equation  explicitly  equivalent  to  Eq,  (23)  was  postulated, 
again  with  2x2  matrices,  and  these  coefficients  were  equated  to  the  coefficients  of  some  two- 
component  subset  of  Bialynicki-Birula’s  four-component  spinor  equation.  The  exact  form  of 
a  free  four-component  Dirac  spinor  in  2D  was  solved  and  used  as  a  guideline  for  choosing 
the  spinor  components.  Here  are  the  positive  and  negative  energy  spinors  in  the  rest  frame 
at  t=0: 

Positive  energy 


1  ^ 

/  0  \ 

- 1 ,  ^  /  E  +  me2 

2  me’ 

0 

j  E  +  me2 

1 

0 

,  u2(p)  =  1 

V  2mc2 

P-C 

E~$~mc2 

o  J 

(32) 
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Negative  energy 


u>3(p)  = 


E  +  me2 
2  me2 


n >  \ 


P+c 


J?+mc2 

l 

o 


*4(p)  = 


E  +  me? 
2  me2 


£+mc2 

0 


0 

1 


It  is  interesting  that  for  infinitesimal  0  values  the  elements  of  the  rows  of  Bialynicki-Birula’s 
D  matrices  have  a  similarity  to  the  large  and  small  magnitude  elements  of  the  fundamen¬ 
tal  orthonormal  Dirac  spinors  in  the  above  shown  ’large  and  small’  representation  and  in 
the  case  of  the  null  lightcone  surface  p^  =  0.  The  resulting  equations  are  inconsistent. 
Thus  to  obtain  the  lattice-gas  representation  of  the  2D  Dirac  equation  using  2x2  matri¬ 
ces  one  has  to  either  guess,  derive  the  explicit  form  of  the  2x2  propagator  solutions  as 
Feynman  [S'c/noe6er,  S.  S.,  1986]  attempted  (and  remains  currently  unsolved),  or  try  to  ’’re¬ 
verse  engineer”  Bialynicki-Birula’s  non-unitary  four  component  spinor  difference  equation. 
Doing  the  latter  did  yield  some  unitary  results.  The  lattice-gas  equation  obtained  is: 


£/$(x,yft)  —  $(%)  2/j  t  t), 

(34) 

where  U  is  a  4x4  matrix  and 

U  =  {SxOiSv  +  GSzdiSy  +  Sx03SyG  +  GSx64SyG)S0 

(35) 

with 

M 

(A1+A2)/2  ~(D1  +  D2)/2\ 
(C1-C2)/2  (B1-B2)/2  ) 

q  _  /  +  A2)/2 

’  U2~\{Cx-C2)I2  ■ 

(Dl  +  D2)/2  \ 
-(B1-B2)/ 2)' 

(36) 

03=  f 

(Bi+ft)/: 2  (C,+C2)/2\ 

-(D,-C2)/2 

a_(-(B>+B2)/ 2 
'  U‘-  UDi-ft)/2 

(C1  +  C2)/  2\ 

(Ai  —  A2)/2  J  ’ 

and, 

Mo  a)- 4-(o  £)•  *- 

and, 

G  =  (o  G)'wlthG=  (l  o) 


/I  0  0  0\ 
0  0  0  1 
0  0  10 
Vo  i  o  o/ 


(37) 


(38) 
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and 


* 


ft 


a  _fCos(me)  iSin(me )  \  a  _(  Oos(me)  iSin(me)  \ 

1  \iSin(me)  Cos(me)  )’  2~\-iSin{me)  -Cos(me) )  ’ 

B  _  f  —iSin(me)  Cos(me)\  _(  -iSin(me)  Cos(me)\ 

1  \—Cos(me)  iSin(me))  ’  2~\Cos(me )  -iSin(me)  J  ’ 

C  —(  iSin{rne)  Cos(me )  \  „  _f-iSin(me)  -Cos(me)\ 

\  Cos(me )  iSin(me)  J  ’  2  \  Cos(me )  iSin(me)  J  ’ 

D  =(Cos(.me )  —iSin(me)\  _(  - Cos(me )  iSin(me)  \ 

1  ^iSm(me)  - Cos(me )  j  ’  2  “  iSin(me )  -Cos(me)  J  ’ 

with  G  a  swap  matrix  for  a  two-component  spinor,  and  and  are  the  usual  streaming 
operators  that  act  on  two-component  spinors.  All  component  matrices  of  U  are  unitary. 
Since  the  streaming  operators  do  not  have  explicit  form  in  this  representation  the  following 
assignments  develop  the  unitarity  of  U. 

The  following  operator  equalities  are  found  to  hold: 

0|04  =  -0x03  (40) 

030|  =  -OxO^. 


By  considering  terms  from  U^U  =  1,  we  are  led  to  determine  the  constants  a  and  b  in: 

0|04  =  -0j03  =  aSyG£fly  (41) 

030i  =  -Oi0\  =  bSxGSl 


Using  Eq.  (41)  in  U  gives  the  following  terms, 


GSX02Sy  =  --S^Sy 

sxo3svg  =  -Hxd3d\o3sy 
Gsxd4syG  =  -\sxd3o\d3sy, 


(42) 
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which  when  substituted  into  U  gives,  with  a=2  and  b=l: 


U  =  -SxOpSySo 
op  =  d3d\d3. 


This  U  is  unitary  and  is  thus  suitable  for  type-II  quantum  computers. 


5.  THE  THREE-DIMENSIONAL  WEYL  AND  DIRAC  EQUATIONS 

The  3D  Weyl  differential  equation  is  obtained  from  Dirac’s  equation,  (17),  with  m=0: 

He 

ihdtip  =  —  a  ■  Vi)  (44) 


There  are  3!  possible  representations,  or  equations,  depending  on  which  components  of  !? 
are  substituted  for  which  components  of  ct.  One  possible  choice  ax  =  az,  ay  =  ax,  and 
Qz  =  <Jy  gives  the  explicit  equation 


dti>i  =  c{-dxipi  -  dytp 2  +  idzfa) 


dtip 2  =  c(+dxip2  -  dy\!)i  -  idzipi). 


(45) 


The  following  unitary  lattice-gas  algorithm  reproduces  this  3D  Weyl  equation  to  first  order 
in  lattice  spacing: 


<f>L (x,y,z,t  +  r)\  _  ( 
4>r(x,  y,z,t  +  T)J  V  3 


o2  sv  a 


q  \  (  4>L(x,y,z,t)\ 
x '  \<i>R(x,y,z,t)  ) 


(46) 


where  Sv,  and  Sz  act  as  the  streaming  operators.  The  C  matrices  are  determined  alge¬ 
braically  by  expanding  the  two  component  spinor  difference  equation  (46)  and  equating  the 
corresponding  matrix  elements  of  equation  (1)  substituted  with  the  2x2  W  matrices  given 
in  Eq.  (12)  of  [Bialynicki-Birula,  /.,  1994],  These  W’s  are  simple  multiples  of  the  2D  W’s 
used  above.  Three  equivalent  matrices  are  then; 


/  —  exp  (ztt/4)  —  exp  (iw /4) 
\^exp(-Mr/4)  —  exp(-irr/4) 


and  Cz  =  C%  =  C\. 


(47) 


Equation  (46)  solves  our  particular  representation,  (45),  of  the  3D  Weyl  equation.  The 
lattice-gas  solutions  to  the  3D  Dirac  equation  may  emerge  from  algebraic  comparison  similar 
to  one  of  those  stated  above.  The  W  matrices  required  are  supplied  in  Bialynicki-Birula  and 
four  component  spinor  lattice-gas  matrices  may  be  required. 


* 
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